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Abstract 

We consider the effect of gold nanospheres of subwavelength sizes on the decay and energy transfer 
rates of quantum systems placed in the proximity of these nanospheres. We find that, for the sphere sizes 
considered in this contribution, the radiative decay rate is barely affected by the presence of the nanosphere, 
whereas the non-radiative decay rate is greatly enhanced due to energy transfer from the quantum system 
to the nanosphere, leading to a strong quenching of the emission of the quantum system. The emission 
wavelength of the quantum emitter and its intrinsic quantum yield play an important role and the impact 
of both has to be considered together when investigating their effect on the non-radiative decay rate. The 
energy transfer process from the emitter to the nanosphere presents a complicated distance dependence, with 
a r~^ regime, characteristic of the Forster energy transfer mechanism, but also exhibiting other distance 
dependence regimes. In the case of a donor-acceptor pair of quantum systems in the presence of a gold 
nanosphere, the donor couples strongly to the nanosphere, acting as an enhanced dipole; the donor-acceptor 
energy transfer rate then follows a Forster trend, with an increased Forster radius. The coupling of the 
acceptor to the nanosphere has a different distance dependence, and it does not follow a Forster-type trend. 
The angular dependence of the energy transfer efficiency between donor and acceptor has a strong dipole- 
dipole trend for small spheres and deviating from it for larger spheres, especially when the donor and 
acceptor are on opposite sides of the sphere. The spectral overlap of the donor emission, acceptor absorption 
and gold nanosphere extinction/scattering shows an interesting trend in that the largest Forster radius is 
obtained when the donor emission and acceptor absorption maxima are somewhat red-shifted from the 
localized surface plasmon peak in the extinction spectrum of the gold nanosphere, being located between it 
and the near-field scattering maximum. 

PACS numbers: 33.80.-b, 42.50.-p, 73.20.Mf 
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I. INTRODUCTION 


The optical properties of quantum emitters (QEs), such as quantum dots and fluorescent dyes, 
near metal nanoparticles (NPs) are a topic of great interest both from the point of view of a fun¬ 
damental understanding of these properties, as well as the technological applications that they 
can lead to, e.g. solar cells,"^^® light-emitting diodes,^subwavelength imaging,^^’^® micro- 
lasers,^^’^^ single photon sources,etc. 

Of particular importance is the influence that the localized surface plasmon (LSP) excited at 
the surface of the metal nanoparticles^^”^^ has on the optical properties of quantum emitters. The 
influence of the LSP on the non-radiative decay rates of quantum emitters placed close to a single 
metal sphere has been extensively studied, both experimentally and theoretically. Different 

models have been proposed for the mechanism of the energy transfer process between the quantum 
emitter and the metal NP, such as the Nanometal Surface Energy Transfer (NSET) model, which 
considers the sphere as an infinite surface and has a distance dependence of the form d~^, and 
the Forster Resonance Energy Transfer (ERET) model, which considers the sphere a polarizable 
dipole, and has a distance dependence of the form Recent theoretical work has shown that 

the distance dependence is more complex, and the NSET and FRET models are only valid in 
certain distance regimes.^^ 

The metal NP is an extended object and therefore has to be modeled in a full quantum elec¬ 
trodynamic formalism. When such a formalism is employed, the distance dependence of the 
non-radiative decay rate of the QE becomes more complex, including r~^, r~^, as well as other 
contributions, in agreement with previous theoretical In addition to the distance 

dependence, we investigate the effect of the emission wavelength and quantum yield, as well as 
the size of the metal NP on the energy transfer rate between the QE and the metal NP. 

Of perhaps even more importance is the process of excitation energy transfer in a donor- 
acceptor pair, one of the main pathways by which energy transfer occurs at the nanoscale. It plays 
important roles in biology,^^^^ nanophotonics (LEDs, nanolasers), microscopy,^^’"^^ sensing,"^^ 
and as optical rulers,^^“^' etc. The process is determined by the dipole-dipole interaction between 
the donor and acceptor and, in the short distance limit, has a r~^ dependence on the donor-acceptor 
separation r, the above-mentioned FRET regime.^^ 

The FRET regime is known to hold for donor-acceptor separations smaller than approximately 
10 nm; for separations larger than this, the energy transfer process is overwhelmed by other de- 
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excitation mechanisms of the donor, such as phonon relaxation or emission into the far-field. In 
a variety of applications, it is desirable to extend the FRET regime to larger donor-acceptor sep¬ 
arations. This would enhance the efficiency of light-harvesting-^^’^'^ and emitting systems,^^’^^ as 
well as increase the range and accuracy of sensing devices"^^ and optical rulers."^^"^^ 

One possible way of enhancing the energy transfer rate in a donor-acceptor pair is by intro¬ 
ducing a “mediator”. The role of the mediator is to increase the probability that energy transfer 
between donor and acceptor will occur. One such mediator is a metal nanoparticle which can sup¬ 
port LSPs, whose main characteristics are both a large confinement and a large enhancement of the 
electromagnetic field. The field of a ESP can thus provide a good mediator for the energy transfer 
process between donor and acceptor, when these are suitably positioned around the nanoparticle. 

The field of the LSP on the metal nanoparticle influences not only the energy transfer rate be¬ 
tween donor and acceptor, but also the other channels which contribute to the total decay rate of 
the donor and introduce new ones, such as absorption of energy by the metal. It is therefore im¬ 
portant to enhance not only the energy transfer rate, but in particular the energy transfer efficiency, 
which is a measure of the competition between the donor-acceptor energy transfer rate and all 
other decay channels of the donor. 

In this contribution we undertake an in-depth investigation of the role of the LSP on small 
spherical Au nanoparticles (Au NPs) in modifying the decay and energy transfer rates and effi¬ 
ciencies of quantum systems placed in the proximity of the Au NP. We begin in Sec. II by laying 
out briefly the theoretical framework used in our investigations. Section III A deals with the decay 
rate of a quantum emitter. In Sec. Ill A1 we validate the model used by simulating experimental 
results obtained in our lab.^^ We then continue with an investigation of the decay rates of a single 
quantum emitter placed near a Au NP and the dependence of these rates on the emission wave¬ 
length and intrinsic quantum yield of the emitter (Sec. Ill A 2), on the distance from the emitter to 
the Au NP (Sec. Ill A3) and on the size of the Au NP (Sec. Ill A 4). In Sec. IIB we consider the 
energy transfer rate between a donor-acceptor pair placed near the Au NP, beginning in Sec. Ill B 1 
with a validation of the model through simulating experimental results obtained in our lab,^^ and 
continuing with an investigation of the dependence of the energy transfer efficiency in the case of a 
single donor-Au NP-acceptor triad on the donor and acceptor distance to the Au NP (Sec. Ill B 2), 
on the relative angular position of the donor and acceptor around the Au NP (Sec. Ill B 3), on 
the spectral overlap between donor emission, acceptor absorption and Au NP extinction spectra 
(Sec. Ill B 4) and on the size of the Au NP (Sec. Ill B 5). Finally, the conclusions are discussed in 
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Sec. IV. 


II. THEORETICAL ERAMEWORK 

Using a Green’s tensor formalism, which we will briefly sketch in what follows, we can calcu¬ 
late the energy transfer rate and efficiency between two point dipoles placed near a metal nanopar¬ 
ticle. 


A. Decay Rates 


Considering a quantum system modeled as a electric point dipole, its decay rate in the presence 
of the Au NP can be related to the power emitted as 


y(r,a)) ^ P{r,co) 

ro{(o) Po{(o) ’ ^ ^ 

where /(r, ®) and 7o(®) are the decay rates of the quantum system in the presence and absence of 
the Au NP, respectively, and P(r, ®) and Po(®) are the power emitted by the dipole in the presence 
and absence of the Au NP, respectively. 

The emitted power in the absence of the Au NP, Po{(o) has the classical expression:^^ 


Po{(o) 


nco'^IJ.Q 


( 2 ) 


The emitted power can be calculated from the integral of the normal component of the Poynting 
vector along a closed surface containing the quantum system: 


P(ro,(y) = ^Re(f)dQ.n E{r,0}) x H*(r, ®). 




(3) 


The total and radiative emitted power can be calculated from this expression by letting the imagi¬ 
nary surface shrink to a point or be infinite. The total power emitted in the presence of the sphere 
is 


P(ro, (o) 



dan-E(r,G)) xH*(r, ®) 


r-5-ro 


E2(ro,6}) 

Eo(ro,®) 



6kc 
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Im(n-6(ro,ro,G))-n) , 


(4) 


4 












where 0(ro, rg, (o) is the Green’s tensor, and the prefactor in the second line is the field enhance¬ 
ment factor, A(ro, (o) = |E 2 (ro, to)/Eo(ro, ®) |^. Similarly, the radiative power emitted in the pres¬ 
ence of the sphere is 


^r(ro,®) 



dQn-E(r,to) xH*(r, to) 


r^oo 


(5) 


These general expressions take, for a spherical NP, specific forms depending on the orientation 
of the transition dipole of the quantum system with respect to the NP. For a radially oriented 
transition dipole, the radiative and total emitted power become 


^r^(ro, to) 
Po{(0) 

P-‘-(ro,to) 

Po(®) 


— A(ro,(o)-^n(n+ l)(2n+ 1) 


Wn (Ps ) (Ps ) 


= A(ro,to) 


l + |l«(n+l)(2«+l)Re(^!^!^ 

^ n \ Ps 


Pj^{ro,(o) _ P-^(ro,tQ) _ Pr^(ro,to) 
Po{(o) ~ Po{(0) Po{(o) 

while the same quantities for a tangentially oriented dipole read 
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Pr (ro,to) 
Po(®) 


= A(ro,to)|j^(2n+l) 

^ n 


Wn{Ps) + f^nCn{Ps) 

2 

1 

V<(Ps)+«nC(Ps) ^ 

Ps 


Ps 


(6a) 

(6b) 

(6c) 


(7a) 


Po(®) 


= A(ro,(o) 


l + ^£(2« + l)Re 

^ n 


bnQiPs) 

Pi 



(7b) 


pJr(ro,to) ^ pll(ro,to) Pj'(ro,to) 

Po{(o) Po{(0) Po{(0) ’ ^ 

In all the above expressions the summation is over the multipole moments n of the sphere, Xj/n 
and represent the first and third kind Ricatti-Bessel functions, Ps = krQ, and a„ and bn are the 
Mie scattering coefficients of the sphere. The relative decay rates can now be calculated from 
Eq. (1). 


B. Energy Transfer Rate 


The energy transfer rate between donor and acceptor quantum systems can be calculated similar 
to the decay rates from the expression 


?bA(rA,rD)PP) _ ^A(tQ) 
7o(®) “ Po{(o) 


( 8 ) 
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where /bA(®) represents the power emitted by the donor D and absorbed by the acceptor A. 

It can be shown that the energy transfer (ET) rate between a donor-acceptor pair modeled as 
point dipoles, is given as^^ 


ybA(rA,rD) = J dm/D(©)7bA(rA,rD, ®) 

0 

oo 

Yd r 

= \^n— / dG)/D((y)|nA-6(rA,rD,a))-nD|^(TA(a)). (9) 

0 

where Td is the donor lifetime in the absence of the acceptor, Fd is the quantum yield of the 
donor, /d(®) is the area-normalized donor emission spectrum, i.e. d(of£){(o) = 1, (7^(0)) is the 
acceptor absorption cross-section and <3 (r ^, re, to) is the Green’s tensor in the particular geometry. 
Finally, nD(A) is a unit vector along the direction of the transition dipole moment of the donor 
(acceptor). 


III. RESULTS AND DISCUSSION 

A. Decay Rates 

A question of some interest and which has also generated some discussion in the literature re¬ 
gards the precise nature of the energy transfer process between a donor quantum system (fluores¬ 
cent dye, quantum dot) and a small Au NP. One candidate mechanism is FRET,^^ which assumes 
that the Au NP can be regarded as a polarizable dipole and the energy transfer process is, therefore, 
that between two point dipoles, with a characteristic r~^ dependence on the donor-acceptor dis¬ 
tance, r. A second mechanism is NSET,^’^^ which views the energy transfer process as occuring 
between a donor dipole and an infinite planar surface, with a characteristic dependence on the 
donor-surface distance, 

Recently, Sukharev et have calculated the distance dependence of the non-radiative de¬ 
cay rates of point dipoles near Ag nanospheres with a diameter of 40 nm. The results of their 
calculations have shown various distance dependence regimes, accommodating both FRET and 
NSET type dependences. In this section we consider the distance, as well as emission wavelength, 
quantum yield, and Au NP size dependence of the energy transfer rate between an emitter and Au 
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(a) Schematic 



(b) Optical properties 



(c) QD-1 


(d) QD-3 


Figure 1: (Color online) Experimental verification of the model used in this paper; (a) schematic of the QD 
monolayer (blue spheres) on top of the Au NP monolayer (gold spheres); (b) normalized emission spectra 
of the two quantum dot samples, as well as the absorption spectrum of the Au NPs. The experimental data 
are taken from Ref. 37 for the quantum dots labeled as (c) QD-1 and (d) QD-3 in that reference. The black 
curve shows our parameter-free simulation. The red dashed curve is a fit of the experimental data with a 
FRET model, while the blue dotted curve corresponds to a fit of the experimental data with a NSET model. 


nanospheres of sizes ranging from subnanometer, when no LSP is supported, to a few hundred 
nanometers, well in the radiative regime. 


7. Experimental Verification 


We first verify our model by considering experimental measurements performed in our lab and 
reported in Ref. 37. Panel la of Fig. 1 shows a schematic of the experimental samples, prepared 
using a Layer-by-Layer technique. A monolayer of Au nanospheres with a radius a = 2.75 nm is 
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separated by a polyelectroyte layer of thickness d (not shown here) from a monolayer of closely 
packed quantum dots (blue spheres). Panel lb presents the optical properties of the quantum 
dots and the Au NPs. The black and red curves represent the emission spectra of two types of 
quantum dots, labeled QD-1 and QD-3 in Ref. 37 and in the rest of this paper, whereas the blue 
curve represents the extinction spectrum of the Au NP monolayer. The photoluminescence (PL) 
spectra and the time-resolved PL decays of these samples were recorded. The quenching of the 
QD emission was calculated from these data as follows. The PL quenching was calculated as 
2 pl = 1 — 7onAu/fQD> where /onAu (^qd) represents the integrated spectral emission of the QD 
monolayer in the presence (absence) of the Au NP monolayer. Similarly, the lifetime quenching 
was calculated as 2 lt = 1 — 't’onAu/'fQD. with TonAu ('fQo) being the PL decay times of the QD 
monolayer in the presence (absence) of the Au NP monolayer. Panels Ic and Id of Fig. 1 show the 
PL and lifetime quenching efficiencies as red squares and blue diamonds, respectively, for the two 
types of quantum dots, QD-1 and QD-3. The agreement between these two sets of experimental 
data, within the error bars, shows that the radiative decay rate of the quantum dots is not modified 
significantly by the presence of the Au NP monolayer. We also show two theoretical fits of the 
experimental data with the FRET model (red dashed curve) and the NSET model (blue dotted 
curve). 

The results of our simulation are shown as a continuous black line in panels Ic and Id. The 
simulations have been performed using the experimental values of the quantum dot properties - 
size, quantum yield, emission spectrum - and Au NP properties - size and concentration in the 
monolayer. As such, the simulation uses no free parameters and one can appreciate the close 
agreement between simulation and experiment. The theoretical quenching efficiency is calculated 
as 


Gth = 


Ynr 

Tr + 'C + rnr’ 


( 10 ) 


where 7r and Ynr are the radiative and non-radiative decay rates of the quantum dots in the presence 
of the Au NP monolayer, while is the intrinsic non-radiative decay rate of the quantum dots, 
when the quantum yield is less than 100%. In this interpretation, Ym is the energy transfer rate from 
the donor quantum dot to the acceptor Au NPs. While both theoretical fits and our simulations 
agree with the experimental data within the error bars, the characteristic distances extracted from 
the FRET fit do not agree with those calculated from the spectral overlap (see 37 for a more in- 
depth discussion). There are many experiments suggesting that the energy transfer mechanism 
from the quantum dot to the Au NP has a more complex distance dependence than these models 
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allow for. We will investigate this distance dependence in section (III A 3). 


2. Emission Wavelength and Quantum Yield Dependence of the Decay Rates 

Now that the theoretical model and simulation procedure have been validated by experimental 
measurements, we consider the simple case of a single emitter placed close to a single Au NR 
We investigate the influence of the emission wavelength of the emitter and its intrinsic quantum 
yield on the quenching efficiency due to energy transfer to the Au NR The emission spectrum of 
the quantum emitter will be modeled as a gaussian distribution with a varying central wavelength 
and a full-width at half-maximum (FWHM) of ~ 45 nm, comparable to the experimental emission 
spectra from Fig. Ic. 

The quenching efficiency of the quantum emitter close to a single Au NP can depend strongly 
on the emission wavelength of the quantum emitter and its intrinsic quantum yield. Fig. 2 shows 
the characteristic 50% quenching distance, Rq, as a function of the emission wavelength of the 
quantum emitter, Aem, and its intrinsic quantum yield, Yq, for several sphere sizes. For the smallest 
Au NP size, a = 0.25 nm (Fig. 2a), the dependence of Rq on Aem and Yq is monotonic, since a Au 
NP of this size does not support a LSP 

When considering a sphere of radius a = 2.5 nm (Fig. 2b), the characteristic distance, Rq, 
attains a maximum value at a wavelength Agm = 527 nm, between the peaks of the extinction and 
the near-field scattering efficiency spectra. The spectral dependence of Rq, for a 100% quantum 
yield is illustrated by the solid black line. Alongside this quantity, we have also included the 
extinction efficiency (red dashed curve) and the near-field scattering efficiency (blue dotted curve). 
The values of these efficiencies are unimportant for the present discussion, and they are shown as 
normalized quantities. 

As the sphere size is increased from a = 2.5 nm (Fig. 2b) to a = 5.0 nm (Fig. 2c), a = 7.5 nm 
(Fig. 2d), a = 10 nm (Fig. 2e) and a = 20.0 nm (Fig. 2f), the dependence of Rq on Aem undergoes 
a blue-shift from Agm = 527 nm (a = 2.5 nm) to Aem = 516 nm (a = 5.0 nm), Agm = 513 nm 
(a = 7.5 nm), Agm = 511 nm (a = 10.0 nm) and Agm = 507 nm {a = 20.0 nm). 

The blue-shift of the maximum of the characteristic distance, Rq, with sphere size is particularly 
visible when comparing panels 2c and 2f. For the smaller spheres, the maximum in Rq lies between 
the maximum of the extinction and the near-field scattering efficiencies (2c). As the sphere size 
is increased, the near-field scattering becomes less important, while extinction (absorption plus 
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(a) a = 0.25 nm 


(b) a = 2.5 nm 


(c) a = 5.0 nm 



400 450 500 550 600 650 700^° 

Emission wavelength, Aem [nm] 



Emission wavelength, Aem [nm] 



Emission wavelength, Aem [nm] 


(d) a = 7.5 nm 


(e) a = 10 nm 


(f) a = 20 nm 


Figure 2: (Color online) Emission wavelength and intrinsic quantum yield dependence of the characteristic 
distance, of a quantum emitter near Au NP of several radii, from (a) a = 0.25 nm to (f) a = 20 nm. The 
excitation wavelength is Aexc = 400 nm. The solid black line corresponds to a fixed value of the quantum 
yield. To = 100%, the dashed red line is the extinction efficiency and the blue dotted line is the near-held 
scattering efficiency for each Au NP. 


far-held scattering) begins to dominate, and the maximum of Rq moves towards the peak of the 
extinction efficiency (2f). 

There are a number of points to be made regarding the joint dependence of the characteristic 
distance, Rq, on the emission wavelength, and the intrinsic quantum yield, Tq. Whereas Rq in¬ 
creases monotonically with the intrinsic quantum yield of the emitter, as expected, its dependence 
on the emission wavelength exhibits a maximum in the vicinity of the LSP peak. It is, therefore, 
important, that in experiments investigating the emission wavelength dependence of the quenching 
efficiency of a quantum emitter, one ascertain whether the intrinsic quantum yield of the quantum 
emitter also changes, as changes in this variable may determine the wavelength of the maximum 
in Rq that one observes. 
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3. Distance Dependence of the Decay Rates 


To investigate the dependence of the non-radiative decay rates on the distance between the 
quantum emitter and the center of the Au NP, r, we write the non-radiative decay rate as: 


ynr(r) 

ro 


r 


( 11 ) 


where Rq is a constant to be determined - the so-called characteristic distance - and n is an expo¬ 
nent, taking the value n = 6 for the FRET model. Taking the logarithm of both sides of Eq. (11) 
above, one can rewrite it as 


log 


ynr(r) 

70 J 


nlogRo — nlogr. 


( 12 ) 


The first derivative with respect to log r of the right-hand side of the above equation will give n as 
a function of the distance from the quantum emitter to the center of the Au NP, n = n{r). 

The emission wavelength of the quantum emitter has been fixed at Agm = 525 nm, with an 
intrinsic quantum yield of To = 100%. 

We have considered the distance dependence of the non-radiative decay rate of a quantum 
emitter in the presence of Au nanospheres with small radii, between 0.25 and 20 nm. Pig. 3 shows 
the exponent n of the distance dependence of the non-radiative decay rate of the quantum emitter 
as a function of the distance of the emitter to the center of the Au NP. Also shown is the energy 
transfer efficiency (dashed lines) in the panels of Fig. 3. As this figure shows, the non-radiative 
decay rate follows a r~^ distance dependence over a range of distances for which the FT efficiency 
is around 50%. This is most clear for the smaller NP, with a radius a = 0.25 nm (Fig. 3a), but it 
is visible for the larger ones as well (Fig. 3b), over a smaller range. Given that the x-axis is 
logarithmic, this range is still quite large. As the radius of the Au NP is increased, however, the 
r~^ distance dependence regime becomes less dominant until it completely disappears for panel 
3f in the figure, for which the radius of the Au NP is 20 nm. One can, therefore, infer a FRET- 
type mechanism for the energy transfer rate from the emitter to the Au NP, when the distance 
of the emitter to the surface of the Au NP is larger than the diameter of the NP. The energy 
transfer rate deviates from a simple dipole pair FRET model when the emitter-Au NP surface 
separation is smaller than the diameter of the Au NP, in which case the higher multipole orders 
in the electromagnetic response of the NP have a non-negligible contribution to the ET rate and 
the NP can no longer be accurately approximated as a point dipole. The characteristic distances, 
Ro, shown in the legends of Fig. 3 are the 50% quenching efficiency distances. The value of 
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(a) a = 0.25 nm (b) a = 2.5 nm (c) a = 5.0 nm 



(d) a = 7.5 nm (e) a = 10.0 nm (f) a = 20.0 nm 

Figure 3: (Color online) Exponent n and the non-radiative decay rate efficiency as a function of distance for 
Au NP with radii (a) 0.25 nm, (b) 2.5 nm, (c) 5 nm, (d) 7.5 nm, (e) 10 nm and (f) 20 nm. 


the quantum yield will not affect the exponent, but it will decrease the Rq distance, according to 


4. Au NP Size Dependence of the Decay Rates 

Finally, we now consider the dependence of the characteristic distance on the size of the Au 
NP. For this investigation, however, we will not consider the distance of the quantum emitter to 
the center of the Au NP, Rq, but to its surface, denoted do- The reason for this change is that, as 
the size of the Au NP is increased, Rq is no longer a good measure of the 50% quenching distance, 
since it is almost equal to the Au NP radius, a. Fig. 4 shows the dependence of do on a for three 
different intrinsic quantum yields of the quantum emitter: (4a) Yq = 100%, (4b) Fo = 10% and 
(4c) Yq = \%. Additionally for each panel, we have also calculated this dependence for several 
emission wavelengths of the quantum emitter, Aem = 400 nm, 450 nm, 500 nm, and 550 nm, 
spanning the LSP range. The log-log scale in Fig. 4 shows that the dependence of do on a follows 
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(a) Yo = 100% 


(b) Yo = 10% 


(c)Yo = l% 


Figure 4: (Color online) Dependence of the characteristic distance, Rq, on the Au sphere radius, a, for 
three different emitter intrinsic quantum yields, (a) Yq = 100%, (b) Yq = 10% and (c) Fq = 1%^ and several 
emission wavelengths of the emitter. Both axes are logarithmic. 


a power law of the form 

do = (Aa)" (13) 

with two regimes. For the first regime, valid at small Au NP radii, the exponent n is in the range 
0.5-0.6, while A is in the range 10-50. For the second regime, valid at much larger Au NP radii, the 
characteristic distance to the surface reaches a saturation regime, where it becomes independent 
of the radius, a. The exponent is then n = 0 and we have A ~ 10 nm (4a), much smaller than 
a. In this latter regime quenching of the quantum system emission by the Au NP is significant 
only extremely close to the Au NP surface. As soon as the quantum system is removed from the 
surface of the Au NP, the emission quenching will no longer be of any importance, as the emission 
properties of the quantum system will be dominated by scattering from the Au NP. 

For the first regime, of small Au NP radii, the dependence of do on a is approximately of the 
form do = \/Aa. As can be seen from the figure, the dependence of do on the emission wavelength 
is rather weak. The dependence on the quantum yield. To follows the dependence mentioned 
in the previous subsection. 


B. Energy Transfer Rates 

Next we consider the energy transfer rates between a donor-acceptor pair placed in close prox¬ 
imity to a Au nanosphere. 
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(a) Schematic 


(b) Optical properties (c) Experimental results 


Figure 5: (Color online) (a) Schematic of the Au nanosphere monolayer (radius 2.75 nm) sandwiched be¬ 
tween the donor QD monolayer (blue spheres) and the acceptor QD monolayer (red spheres); (b) Emission 
spectrum of donor QD monolayer (black solid line), absorption spectrum of the acceptor QD monolayer 
(red dashed line) and extinction of Au NP monolayer (blue dotted line) (c) Energy transfer efficiency versus 
the center-to-center distance for monolayers of quantum dots separated by a monolayer of Au nanospheres 
of radius 2.75 nm. The donor-Au NP separation is kept fixed at 3 nm. See reference 57. The black solid 
line represents our simulation results. 


1. Experimental Verification 


We validated the model by comparing it with experimental measurements performed in our 
laboratory on energy transfer rates between monolayers of donor and acceptor quantum dots, sep¬ 
arated by a monolayer of Au nanospheres.^^ The schematic of the structure is presented in Fig. 5a, 
where the donor quantum dots are represented as blue spheres in the bottom monolayer, the ac¬ 
ceptor quantum dots are represented as red spheres in the top monolayer, and the Au NPs are 
represented as gold spheres in the middle monolayer. The separation of the different monolayers 
is controlled by the use of polyelectrolyte layers of known thickness (not shown in the schematic). 
Fig. 5b shows the emission spectrum of the donor quantum dot monolayer (black continuous line), 
the absorption spectrum of the acceptor quantum dot monolayer (red dashed line) and the extinc¬ 
tion spectrum of the Au NP monolayer, showing the LSP peak (blue dotted line). This panel shows 
a good overlap between the different spectra, an essential feature for optimizing the energy transfer 
process. Finally, fig. 5c shows the calculated energy transfer efficiency in such a system, together 
with the experimental measurements. In this multilayer structure, the donor-Au NP separation is 
kept fixed at 3 nm, and the acceptor-Au NP separation is varied. The simulation uses experimental 
values of all parameters and is, therefore, parameter-free. The experimental results are shown as 
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red discs, together with the associated error bars, while the results of the simulations are presented 
as the black solid line. A good agreement between the model and experiment is observed. We feel 
justified, therefore, to now use this model in the following investigation of the distance, angular 
position, spectral overlap and Au NP size dependencies of the energy transfer in the presence of 
Au nanospheres. 


2. Distance Dependence of the Energy Transfer Efficiencies 


Starting with this section, we consider a single donor-Au NP-acceptor triad and investigate the 
donor-acceptor energy transfer. In this, as well as the following sections, the intrinsic quantum 
yield of the donor is, Yq = 100%, the donor emission spectrum is modelled as a gaussian with a 
FWHM ~ 45 nm, while we use the experimental absorption spectrum of the acceptor shown in 
panel 5b. 

We first investigate the dependence of the energy transfer efficiency between a donor-acceptor 
dipole pair on their separation. The donor and acceptor are placed on opposite sides of a Au NP of 
various radii, forming a linear triad. Fig. 6 shows the distance dependence of the energy transfer 
efficiency from donor to acceptor for the case when the distance from the donor (6a, 6c, 6e) or 
the acceptor (6b, 6d, 6f) to the Au NP is constant, while the position of the other member of 
the pair is varied. Each panel of Fig. 6 shows the energy transfer efficiency as a function of the 
donor-acceptor center-to-center separation, tda. for several Au NP radii, a = 2.5 nm for 6a and 
6b, a = 5.0 nm for 6c and 6d, and a = 10 nm for 6e and 6f. Panel 6a shows, in the main plot, 
the energy transfer efficiency between the donor and acceptor in the absence of the Au NP (black 
circles), when the donor is dX dD = 2> nm from the surface of the Au NP (red squares) and when 
the donor is at Jd = 6 nm from the surface of the Au NP (blue diamonds). Additionally, fits of the 
calculated data with a function of the form 


^7(^a) 


1 



(14) 


are shown as lines for each data set. The data sets have been fitted to the curves in the region 
where r\ (^da) = 50% and then the fitting curves have been extended over the entire data range. 
The legend shows the fitting parameters Rq and n. The inset in the panel shows the same data sets 
on a log-log scale, together with a different set of fitting curves. The procedure used to obtain this 
second set of curves has been to fix the value of the n parameter to n = 6 (corresponding to FRET) 
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(a) a = 2.5 nm, fixed donor distance (b) a = 2.5 nm, fixed acceptor distance 



Center-to-center separation, Tda [nm] 

(c) a = 5.0 nm, fixed donor distance 



Center-to-center separation, roA [nm] 


(d) a = 5.0 nm, fixed acceptor distance 




(e) a = 10.0 nm, fixed donor distance (0 ci= 10.0 nm, fixed acceptor distance 

Figure 6: (Color online) Energy transfer efficiency between a donor-acceptor pair on opposite sides of the 
Au NP, as a function of donor-acceptor separation, tda, when the donor-Au NP distance is kept fixed (a), 
(c), (b) and the acceptor-Au NP distance is kept fixed (b), (d), (f); The Au NP radius is a = 2.5 nm in (a) 
and (b), a = 5.0 nm in (c) and (d), and a= 10 nm in (e) and (f). The insets show the same data on a log-log 
scale. 
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and fit the data sets in the region where this power law holds, i.e. at relatively large donor-acceptor 
separations. The legend in the inset shows the Rq parameter extracted from these fits. Panel 6b is 
analogous to 6a, except that now the distance from the acceptor to the surface of the Au NP is kept 
fixed at = 3 nm (red squares) and Ja = 6 nm (blue diamonds). Panels 6c-6f are similar, except 
for different radii of the Au NP. 

Beginning with panel 6a, for which the donor-Au NP distance is kept fixed, it is evident that the 
energy transfer process between donor and acceptor follows a dependence on their separation, 
and is, hence, a FRET process. It is, however, a modified FRET process, since the characteristic 
distance or Forster radius, Rq, takes on values dependent on the donor-Au NP distance and larger 
than the Forster radius in the absence of the Au NP, which is Rq = 10.0 nm. Thus, when the 
donor-Au NP distance is = 6 nm, the Forster radius increases to Rq = 181 nm and at = 3 
nm, it is Rq = 28.8 nm. 

When considering panel 6b, for which the acceptor-Au NP distance is kept fixed, the situation 
is quite different. The main plot shows an increased characteristic distance Rq = 14.8 nm at «iA = 3 
nm and Rq = 17.3 nm at Ja = 6 nm, but the values of n one obtains, n = 17.5 at <iA = 3 nm and 
n= 13.6at<7A = 6 nm, are very different from n = 6, making it clear that this is no longer a FRET 
process, and the reduction of the energy transfer efficiency with distance is much more dramatic. 
We see in the inset that, were one to insist on a fit with a FRET model, that fit would be valid 
only for donor-acceptor separations larger than 20 nm. The Rq extracted from this fit differs only 
negligibly from the free-space Rq and is reduced, rather than enhanced. 

The lessons to be drawn from the first two panels of Fig. 6 are interesting: a fixed donor-Au 
NP separation leads to a FRET-type behavior of the energy transfer efficiency, with an increased 
characteristic distance Rq, which can be almost tripled. This suggests that such a system could 
be profitably used as a spectroscopic ruler: the FRET-type behavior assures a consistent distance 
dependence over large ranges, while the increased characteristic distance Rq extends the functional 
range of the ruler. A fixed acceptor-Au NP system, on the other hand, is not suitable as a spectro¬ 
scopic ruler because, even though the characteristic distance Rq is increased, there is no consistent 
distance dependence of the energy transfer efficiency. 

Considering now the second (6c and 6d) and third (6e and 6f) sets of panels in Fig. 6, the first 
things to notice is an increase in the characteristic distance, Rq, in the legends of the main panels, 
and a progressively more pronounced deviation of n from the FRET value of n = 6. These are a 
result of increasing the radius of the Au NP from a = 2.5 nm to a = 5.0 nm, to a = 10.0 nm, with 
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the effect that the range of donor-acceptor interaction is increased, and the higher-order multipoles 
of the Au NP contribute to this interaction. Secondly, the characteristic distances in the insets of 
each panel, extracted from a pure FRET fit, present an increase with the sphere size for a fixed 
donor-Au NP separation, but not for a fixed acceptor-Au NP separation. This happens because the 
FRET fit is valid only at large overall donor-acceptor separations, where the ET efficiency reverts 
to its value in the absence of the Au NP. 

3. Angular Position Dependence of the Energy Transfer Efficiencies 

In the previous subsection, the donor and acceptor were placed diametrically opposed around 
the sphere, i.e. at antipodes, and the energy transfer efficiency was investigated as a function of 
their distance to the surface of the sphere alone. We now investigate the behavior of the energy 
transfer efficiency when the angular position of the donor and acceptor around the sphere is varied. 
Fig. 7 shows polar plots of the ET efficiency for a donor-acceptor pair close to a Au NP of radius 
a = 2.5 nm (left column), a = 5.0 nm (middle column) and a = 10 nm (right column). The position 
of the donor is fixed on the x-axis at <7 d = 3 nm (top row), d\y = 6 nm (middle row) and d\y = 9 
nm (bottom row) from the surface of the Au NP, while several acceptor distances are considered 
(see legend in Fig. 7). The emission wavelength of the donor is Agm = 525 nm, and an average is 
performed over donor and acceptor dipole orientations. The thin curves in each panel represent 
the ET efficiency in the absence of the Au NP, i.e. in free-space, for exactly the same parameters. 

For the first row of Fig. 7, comprising panels 7a, 7b, and 7c, the distance of the donor from 
the surface of the Au NP of different sizes is <7 d = 3 nm. In this case, the ET efficiency around 
the Au NP is very close to 100%, for all three sphere sizes, and all acceptor distances to the 
surface of the Au NP, though it can be seen that the ET efficiency begins to decrease for the largest 
acceptor distance, dp^ = 9 nm. This effect is more pronounced for the larger Au NP, for which the 
donor-acceptor separation is larger. 

As the donor distance to the surface of the Au NP is increased, the angular dependence of the 
ET efficiency is different for Au NPs of various sizes. In the case of the smallest Au NP with a 
radius a = 2.5 nm (left column in Fig. 7), the angular dependence of the ET efficiency resembles 
the ET efficiency in the absence of the Au NP, albeit with larger overall values. Furthermore, 
the results obtained with this full method are practically indistinguisable from those obtained by 
treating the Au NP as a point dipole (data not shown). This strongly suggests that the FRET-type 
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(a) a = 2.5 nm, Jd = 3 nm 


(b) a = 5.0 nm, Jd = 3 nm 


(c) ( 2=10 nm, Jd = 3 nm 


90° 


90° 


90° 



(d) a = 2.5 nm, Jd = 6 nm 


(e) a = 5.0 nm, Jd = 6 nm 


(f) a= 10 nm, Jd = 6 nm 



Figure 7: (Color online) Polar plot of the ET efficiency between a donor located at Jd = 3 nm (top row), 
dj) — 6 nm (middle row) and Jd = 9 nm (bottom row) from the surface of the Au NP on the x-axis and 
an acceptor located at different distances, Ja, from the surface of a sphere of radius a = 2.5 nm (left 
column), a — 5 nm (middle column) and a = 10 nm (right column). The emission wavelength of the 
donor is Aem = 525 nm. The orientations of the donor and acceptor dipoles are isotropic. The thin lines 
correspond to the same arrangements of donor and acceptor as the thick lines, except that they are calculated 
in the absence of the Au NP 
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(a)a = 2.5nm (b)a = 5.0nm (c)a = 7.5nm 

Figure 8: (Color online) ET rate enhancement as a function of the central wavelength of the donor emission 
spectrum for several positions of the donor and acceptor and three sizes of the Au NP: (a) a = 2.5 nm, (b) 
a = 5 nm and (c) a = 7.5 nm. 


interaction we have seen being valid in the previous section on the distance dependence of the ET 
rate, also describes the angular dependence in the case of the smallest sphere, a = 2.5 nm. 

For Au NP of larger sizes (the middle, a = 5 nm, and right columns, a = 10 nm, in Fig. 7), 
the angular behavior of the ET efficiency begins to deviate from a FRET-type. This is especially 
notable in the appearance of a bulge at an angle 7t, which contrasts with the dip observed in 
the free-space ET efficiency at the same angle. Moreover, when treating the Au NP as a point 
dipole, this bulge is absent, strongly suggesting that it is due to higher order multipoles of the Au 
NP. A final piece of evidence for this hypothesis is the fact that the bulge is less pronounced or 
dissappears when either the donor or acceptor is moved away from the surface of the Au NP and 
can no longer couple to these higher order modes. 


4. Spectral Overlap Dependence of the Energy Transfer Rates 

We now consider the dependence of the energy transfer rate and efficiency on the spectral 
overlap between the donor emission, acceptor absorption and the ESP peak of the Au NP. To this 
end we fix the position of the donor and acceptor at several distances from the surface of the Au 
NP, on opposite sides. We then consider gaussian distributions with a fixed width of 45 nm for 
the donor emission spectrum and we calculate the energy transfer rate as a function of the central 
wavelengths of these distributions. The results are shown in Fig. 8 for three sizes of the Au NPs, 
a = 5,10,15 nm. As the figure shows, the ET rate enhancement peaks between 500 nm and 550 
nm for all sizes. For the smallest size, a = 2.5 nm (panel 8a), the ET enhancement peak is around 
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530 nm when both donor and acceptor are very close to the An NP (3 nm) and it blue-shifts to 
520 nm when the donor and acceptor are moved away from the Au NP. As the size of the Au 
NP is increased to a = 5 nm (Fig. 8b) and a = 7.5 nm (Fig. 8c), the FT rate enhancement peak 
is around 525 nm and is not noticeably dependent on the position of either the donor or acceptor. 
For the case of the smallest Au NP, Fig. 8a with a = 2.5 nm, the LSP is around A = 500 nm. The 
maximum FT rate enhancement for this case, however, occurs between the LSP peak at 500 nm 
and the near-field scattering peak at 540 nm (see Fig. 2 for the spectral shapes of the extinction and 
the near-field scattering efficiencies of Au NPs of various sizes). This reflects the influence of the 
near-fleld of the Au NP, which has evanescent components not visible in far-field measurements, 
but which contribute to FT close to the Au NP. 

Another feature of the plots in Fig. 8 is the broadening of the spectra with decreasing Au NP 
size. This is clearly seen in Fig. 8 and is a consequence of the broadening of the LSP peak due to 
a surface scattering contribution in the dielectric function of Au spheres of sizes below 50 nm.^^ 

5. Au NP Size Dependence of the Energy Transfer Efficiencies 

Finally, we investigate the dependence of the FT efficiency on the size of the Au NP. Fig. 9 
shows the sphere radius dependence of the FT efficiency for several distances of the donor from 
the surface of the Au NP, dy, = 1 nm (panel 9a), dy, = 3 nm (panel 9b), dy, = 5 nm (panel 9c) 
and d\y = l nm (panel 9d), and several distances of the acceptor from the surface of the Au NP 
(<7 a = 1 nm to <7 a = 9 nm). The donor and acceptor are on opposite sides of the Au NP. When 
the donor distance to the surface of the Au NP is small, as in panel 9a, the FT efficiency is almost 
100%, for sphere radii up to 100 nm, when the acceptor distance to the surface is also relatively 
small. As the acceptor distance to the surface is increased, panel 9a shows the development of 
a peak in the FT efficiency around and above 100 nm. This extra peak is not present when the 
Au NP is treated in the dipole approximation; indeed, only when one takes into account multipole 
orders above the tenth, does the peak appear. It is therefore reasonable to attribute this peak to the 
higher order multipoles of the Au nanosphere. This is corroborated by the fact that the height of 
the peak decreases rapidly as the distance of either the donor or acceptor to the surface of the Au 
NP is increased. This suggests that the effect belongs to the near-field of the Au NP, which is in 
accordance with the effect being due to the higher multipoles of the Au nanosphere. Increasing 
the donor distance to the surface of the Au NP, as in panels 9b to 9d, the behaviour of the FT 


21 



0.1 1 10 100 1,000 0.1 1 10 100 1,000 


Sphere Radius, a [nm] Sphere Radius, a [nm] 

(a) r/o = 1 nm (b) r/o = 3 nm 



Figure 9: (Color online) ET Efficiency as a function of the sphere radius, a, for several acceptor positions 
and four donor positions, (a) dD = I nm, (b) do = 3 nm, (c) — S nm and (d) d^ = 7 nm. 


efficiency changes quite dramatically, especially for very small sphere radii. When the radius of 
the Au NP is significantly smaller than the distance of either the donor or acceptor to the surface 
of the Au NP, the ET efficiency can be considerably reduced. Our simulations show, in fact, that 
in the limit a^O, the ET efficiency in the presence of the Au NP is equal to the ET efficiency in 
free-space (data not shown). 


IV. SUMMARY AND CONCLUSIONS 

In this contribution we have investigated the quenching of the emission of a quantum emit¬ 
ter placed close to Au NPs. We have employed a Green’s tensor formalism, which we have first 
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validated by simulating experimental data obtained in our lab, for monolayers of QDs on top of 
monolayers of Au NPs. We have found good agreement between this parameter-free simulation 
and the experimental measurements. Subsequently, we have employed the Green’s tensor formal¬ 
ism to investigate the dependence of the quenching efficiency of a single QE close to a single Au 
NP on several parameters, such as the emission wavelength, intrinsic quantum yield of the QE, the 
size of the Au NP and the QE-Au NP distance. 

We have considered the dependence of the quenching efficiency on the emission wavelength 
and intrinsic quantum yield of the QE together, since these two parameters are often modified 
in tandem in experiments. We have found that the characteristic distance of the QD-Au NP pair 
increases monotonically with the quantum yield of the emitter, whereas the dependence on the 
emission wavelength of the emitter has a rather pronounced peak located between the ESP wave¬ 
length of the Au NP and the near-field scattering efficiency maximum. When considering the 
dependence of the quenching efficiency - or characteristic distance - experimentally, it is there¬ 
fore important to take into account both of these dependences, as not to do so can give artifacts, 
such as an apparent shift of the maximum quenching efficiency with the emission wavelength of 
the emitter. 

We have also investigated the distance dependence of the quenching efficiency, mainly to as¬ 
certain under which conditions the ERET model for the QE - Au NP interaction applies. Our 
simulations show that the FRET model provides a good approximation of the interaction, espe¬ 
cially close to the ESP wavelength and for NP sizes below 5 nm. For larger spheres, not only does 
the exponent n deviate from the value of 6 associated with the FRET model, but the fit becomes 
poorer, suggesting that the single power law used in the fit no longer accurately applies. 

We must mention here that our model does not reproduce the experimental results obtained 
by the Strousse group'“^ for the energy transfer from a quantum emitter, a fluorescent dye in this 
case, to a Au NP with a diameter of approximately D = 1.5 nm. These authors have found that 
the distance dependence of the quenching efficiency is better described by the NSET model, with 
a d~^ distance dependence. Evidently, more research is required to elucidate this discrepancy. 

Finally, considering the behavior of the quenching efficiency as a function of the size of the Au 
NP, we have seen that the characteristic distance to the Au NP surface is determined by the extinc¬ 
tion cross-section of the NP. Thus, when the size of the Au NP is much smaller than the emission 
wavelength of the QE, the extinction cross-section is dominated by absorption and the character¬ 
istic distance can be 2-3 times larger than the NP radius, following a power law dependence on 
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the Au NP radius. In the size regime where scattering dominates the extinction cross-section, this 
distance becomes constant and independent on the Au NP radius. Quenching, in this latter case, 
occurs mainly very close to the NP surface and the QE emission can be expected to be enhanced 
and not quenched further away from the surface. 

We have also investigated the energy transfer between a donor and acceptor. The energy transfer 
efficiency is much more sensitive to the donor-Au NP distance, than to the acceptor-Au NP dis¬ 
tance. Furthermore, the distance dependence of the energy transfer efficiency for fixed donor-Au 
NP separations suggests that the donor-acceptor interaction mediated by the Au NP is essentially of 
Forster type, following a r~^ dependence on the donor-acceptor separation r over a large distance 
range, but with an increased characteristic distance. This is not the case when the acceptor-Au NP 
separation is fixed. In this case the energy-transfer efficiency does not follow a Forster dependence 
and it is influenced by the presence of the Au NP only for small donor-Au NP separations. 

The positioning of the donor and acceptor around the sphere can also influence the energy 
transfer, especially for smaller Au NPs and donors and acceptors placed very close to the Au 
NP. For small Au NP sizes, the angular distribution of the FT efficiency is similar to that for a 
dipole pair, albeit with an increased characteristic distance. As the size of the Au NP is increased, 
however, there is a marked deviation from the dipole pair behavior at an angle of K. This can be 
attributed to the contribution of the higher order multipoles of the Au NP to the energy transfer 
process. 

When investigating the dependence of the FT transfer efficiency between a donor-acceptor pair 
placed near the Au NP on the emission wavelength of the donor, we have shown that the largest FT 
rate does not occur at the LSP peak, but is red-shifted close to the peak in the near-field scattering 
efficiency. This can be attributed to the fact that, when the donor-acceptor pair is placed in the 
near-field of the Au NP, the evanescent components of the electric field can act as a strong channel 
for the FT process. This is more clearly apparent when one considers the FT rate enhancement 
dependence on the donor-acceptor distance to the Au NP. In this case a gradual blue-shift towards 
the LSP wavelength occurs when the distance is increased and, hence, the evanescent component 
contribution becomes negligible. 

Finally, an interesting behaviour of the FT efficiency as a function of sphere radius has been 
obtained from our simulations. For relatively small distances of the donor and acceptor to the 
surface of the Au NP, in addition to the large FT efficiency obtained for relatively small radii, 
there is an extra peak in the FT efficiency when the radius is around 100 nm. We have shown that 
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this peak is a consequence of the higher multipole orders of the Au nanosphere contributing to the 
energy transfer process, for large enough spheres. 

In this contribution we have undertaken a thorough investigation of the effects of a Au NP 
of subwavelength sizes on the decay and energy transfer processes of quantum systems placed 
in its vicinity. The increased characteristic distance provided by the Au NP can have interesting 
applications in improved sensing and as an extended spectroscopic ruler. 
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